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We analyze the spin-orbit terms in multi-subband quasi-two-dimensional electron systems, and 
how they descend from the bulk Hamiltonian of the conduction band. Measurements of spin-orbit 
terms in one subband alone are shown to give incomplete information on the spin-orbit Hamilto- 
nian of the system. They should be complemented by measurements of inter-subband spin-orbit 
matrix elements. Tuning electron energy levels with a quantizing magnetic field is proposed as an 
experimental approach to this problem. 
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Spin-dependent phenomena in semiconductors was one of the favorite research themes for Vladimir Idelevich Perel' 
since the early 1970s.— From the mid-1980s his interests shifted towards spin-related effects in low-dimensional systems, 
starting with optical orientation and polarization properties of hot photoluminescence in quantum- well structures, 
which was closely connected to experiments carried out by the group of D.N. Mirlin.^ Some of this work was done 
together with one of us (MEP).— In the series of more recent papers, the transition from the two-dimensional to 
quasi-three dimensional case was considered^ Several latest publications of Vladimir Idelevich were focused on spin- 
dependent tunneling and the role played in it by spin-orbit interactioni^ This has defined the subject choice for our 
contribution to the special issue devoted to his memory. 

The spin-orbit interaction in semiconductors has been widely discussed recently in relation to some proposals of 
spin-electronic and quantum-computing devices. It is of considerable physical interest in itself, because due to strong 
gradients of atomic potentials within the crystal unit cell the spin-orbit terms in the effective-mass Hamiltonian are 
often greatly enhanced with respect to those of a free electron.^ In addition, the reduced crystal symmetry brings 
about new spin-orbit terms unknown for free particles, the so-called Dresselhaus terms.— 

In two-dimensional systems, spin-orbit effects are known to be even stronger than in bulk semiconductors. In 
particular, in the effective-mass two-dimensional (2D) Hamiltonian there appear spin-orbit terms which are linear 
in the 2D wave vector k. They may exist in structures where the spatial inversion symmetry is broken. There are 
so-called bulk inversion asymmetry (BIA) terms, which appear on averaging the bulk Dresselhaus terms over the 
envelope function of the corresponding size-quantization level, and structure inversion asymmetry (SIA), or Rashba, 
terms. S The latter are believed to exist in asymmetric quantum wells (QWs); there is plenty of experimental evidence 
of their existence in specific structures, but apparently no agreement has been reached as to how they descend from 
the bulk spin-orbit Hamiltonian. Some authors argue that they are entirely due to interfacial effectsi^ The most 
consistent theoretical treatment of this problem was carried out by Gerchikov and Subashiev^'^ (see also a more recent 
paper by Winkle»ii). However, even these papers do not give explicit answers to questions arising when one attempts 
to devise experiments aimed at determination of spin-orbit parameters or to engineer structures with controllable 
spin-orbit effects. 

In our opinion, for a full understanding of the spin-orbit effects in nanostructures it is necessary to take into 
consideration inter-subband spin-orbit coupling in multi-subband quantum-dimensional structures, which present a 
d ' natural bridge between 3D and 2D semiconductor structures. In this paper we analyze the effect of the spin-orbit 
interaction on the energy spectrum of many-subband QWs, and show that applying a quantizing magnetic fields 
to such a system presents a way for experimental determination of all the relevant parameters of the spin-orbit 
interaction. Implications for optical and phonon spectroscopy are discussed. 

We start with the phenomenological expression for the spin-orbit Hamiltonian of the conduction band of a com- 
positionally homogeneous bulk semiconductor in the envelope-function approximation. It has the following general 
formi^ 

Hso = (hs/A • S)+ay([k X VV] ■ S), (1) 
where the "Dresselhaus field" hsiA^ existing in non-centrosymmetric crystals, is a pseudovector that is an odd function 



2 



of the components of the electron wave vector k, and V is the electrostatic potential energy. The second term in 
Eq. ID) has the same form as the spin-orbit Hamiltonian of a free electron; however, the constant ay is not equal to the 
vacuum spin-orbit constant ayac — (4TOgC^), where mg is the bare electron mass. As electrons in semiconductor 
crystals are actually subjected to strong potential gradients within the crystal unit cell, the spin-orbit interaction is 
enhanced by a factor of approximately rrieC^ / Eg, where Eg is the semiconductor band gap^i^ This huge enhancement 
of the spin-orbit interaction has allowed, in particular, observation of spin-dependent currents due to anisotropic 
scattering of electrons by impurity centers (a solid-state analog of the Mott effect) in GaAsJ^ii^iii^ 

In nanostructured semiconductors based on solid solutions like Ga^jAli-jjAs, not only the electrostatic potential, 
but also the composition x and, respectively, all the parameters of the band structure, may depend on the coordinates. 
One can therefore introduce a new phenomenological spin-orbit term proportional to Vx. Since all the band energies 
are, to the first approximation, linear in x, we shall write this term for electrons as ax([k x VEq] ■ S). Here, VEc is 
the "variable-gap" field that affects charge carriers in structures with gradients of composition.^^ 

Using the standard 8x8 kp method of calculation of spin-orbit splitting, analoguos to that used in Refs. [ol lTGlfTTllTsl . 
Gerchikov and Subashie^ii^ have shown that the spin-orbit term in the conduction band Hamiltonian can be expressed 
through a gradient of an effective "spin-orbit potential" xc- 

Hso = (hB7A+[k X Vxc]) • S, (2) 

where xc is a function of the energy positions of the extrema of the conduction band (Fg) and the uppermost valence 
bands (Fg and Tj): 

3[E-Ev{x)][E-Evix) + A{x)y ^' 

Here, P is the momentum matrix element between S and P Bloch states, E is the electron energy, Ey (x) is the energy 
position of the top of the valence band, and A(x) is the spin-orbit energy splitting. As noted in Ref.llCi Eq. ^ is true 
even in interface regions, where Vxc can be expressed in terms of delta functions. In the following, we shall consider 
only structures where the energy of size quantization is much less than the bandgap Eg . This condition assumes that 
variations of both V and Ec are much less than Eg over the region where size-quantized wave functions are mainly 
concentrated. This class of structures includes wide quantum wells or heterojunctions with high barriers, but in this 
case interfacial regions, where the electron probability density is small but the gradient of Ec is very large, should 
be considered separately (see discussion below). For such structures, one can easily obtain from Eq (2) of Ref. [13 
(compare also Eqs. (7) and (10) of Ref. [Tsh the following expressions for ay and ax- 
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ax 
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where m = 3h'^P^Eg{Eg + A)/2{Eg + 2 A) is the effective mass of the conduction-band electron. For example, using 
band parameters and their dependence on composition in Gai_2.Ala;As from Ref. tl9|, one obtains ay = 4.1 x 10~^^ 
cm~^, and ax = —3.1 x 10~^^ cm~^ near a; = 0. For the class of structures we consider, spin-orbit constants, as well 
as the effective mass, can be approximately treated as spatially invariable parameters. 

We can now analyze how bulk spin-orbit terms transform as a result of size quantization in QWs. Due to trans- 
lational invariance in the QW plane, the general form of the matrix element of the spin-orbit Hamiltonian between 
electron eigenfunctions in the QW is: 



mki/i 



H 



so 



n\.2v) = (h™„(k) • S)^^5(k - ki)J(k - ka). 



(5) 



where m and n enumerate size-quantization subbands, ki and k2 are electron wave vectors, and /i and v are spin 
indices. 

Let us first consider the SIA (Rashba) terms. If all the potential gradients are normal to the QW plane (we denote 
the corresponding unit vector as n), the spin-orbit field h„i„(k) takes the form: 



'^mn* C*^) — -^Vin [n X k] 



(6) 
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where 

+00 



-00 

+ 00 



= ay y 'fmiz)'fn{z)^dz + ax J *™(z)*„(z)^dz + ^ /^„. (7) 

—00 —00 ^ 

Here, vE'm(z) and ^n(z) are envelope functions of the mth and nth subband in the z-direction, respectively. The values 
of av and ax correspond to the bottom of the QW. The last term is an interfacial contribution, which arises because 
at interfaces sharp gradients of Ec coincide with abrupt changes of xc] summation is taken over all the interfaces. 
Since wave functions are continuous at interfaces, the overall contribution of interface regions to the matrix element 
^nrn^ can be written as '^m{zi)'i'n{zi){xc{zi+) — xc{zi~)), where Zi+ and Zi_ denote positions immediately to the 
right and to the left of the interface, respectively. The term is therefore equal to the difference between this 
expression and the contribution of the interface to the integral in the second term of Eq. ([7]) : 

I'mn = '^m{zi)^„{zi) [xc{zi+) - Xc{Zi-)] - (z^ )^'„ (z^ )ajf [Ec{zi+) - Ec{Zi^)] = 

= ^'mizi)'^n{zi){ai^ ax)[Ec{zi+) - Ec{zi^)], (8) 

where di — [xc{zi+) — Xcizi-)] / [EciziJ^) — Ec{zi-)] has the meaning of an effective constant ax, renormalized 
due to abrupt changes of the band structure at the interface. For example, for the GaAs/Ga^ gAlo.4As interface, 
di = —2.5 X 10^^^ cm-^. Since the bulk constant ax = —3.1 x lO"""^^ cm^, the difference \di — ax\ is much smaller 
than \ax\, and we can conclude that the interfacial corrections are not significant in the GaAs/Gaj^ ,j,Ala;As system. 

As the QW localization potential for electrons is U{z) — V{z) + Ec{z), Eq. ([7]) can be rewritten in the following 
way: 

H-oo +00 

'i>,n{z)^n{z)—dz + {av - ax) / *,„(z)*„(z)— dz + ^ /;,„. (9) 

— 00 —00 * 

Making use of the fact that ^ — ~Pz, where pz is the z-component of the electron momentum, we obtain*^ 

+00 +00 +00 

^^{z)-i'„{z)—dz^-UEn-E,n) / •^,n{z)pz-^n{z)dz ^ {E^- E„) / (z) — (z)dz, (10) 

dz h J J dz 

— oo — oo 

where _E„ and Em are energy levels corresponding to eigenfunctions 'l'n(z) and 5'm(z). 

Considering the second term in Eq. ([9]), we recall that £ = —{\/e)/{dV/dz) is an electric field which is constant 
across the quantum well unless there are electric charges inside the well. Modern nanostructure technology usually 
avoids placing impurity atoms inside the quantum well, so an inhomogeneity of £ may normally arise only due to 

screening by the electron gas in QWs containing free electrons. If the concentration of two-dimensional electrons is 

+00 

small, the electric field inhomogeneity can be neglected. For a constant £, J '$rn{z)'^n{z)£dz ~ £S„in. In this case, 

— C30 

we come to the following expression for A^r^- 



/d 
*™(z)— «'„(z)dz - {av - ax)e£Smn + I^n- (H) 

-00 * 
The value A^^"^ = —{av — ax)e£ + ^Inn coefficient of the Rashba term in the nth subband. One can see 

i 

that, contrary to the widespread opinion, the Rashba term is not entirely due to interfacial effects; in fact, simple 
estimations show that the interfacial term is much smaller than the bulk contribution proportional to the electric 
field and to the difference of the spin-orbit constants av and ax- A parabolic quantum well formed by modulation 
of composition is a good example illustrating this property of SIA spin-orbit terms. There are no interfaces in the 
parabolic quantum well. Moreover, applying an electric field in this case does not change the shape of the localizing 
potential, which remains symmetric: 

22 9 22 

, , rabj z „ muj , ,2 muj zf, 

U{z)^— £ez=^(z-zo)' (12) 
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where zq = EejmuP' . However, according to Eq. there exist both intra- and intersubband SIA terms: 



,s/A _ ASIA imhw^n 



^n-l,n — ^n,n-l — "X 1 

A^'n = -{av-ax)e£. (13) 

Intra-subband SIA terms are zero in the absence of electric fields, while inter-subband SIA matrix elements are always 
present and do not depend on the electric field. Remarkably, bandgap gradients do not contribute to the Rashba 
term: one needs to apply an electric field to the structure to produce it. These gradients are needed, however, to 
confine the electron in the z-direction; if the electron was confined by a purely electrostatic potential, Rashba terms 
would also be absent. This fact has been noted by Gerchikov and Subashiev.^'^ 

Another instructive example is a symmetric rectangular quantum well in an electric field. Since in this case all the 
changes of composition x are concentrated at interfaces, the interfacial effects can be exactly accounted for by replacing 
ax with di in Eq. ()lip : this substitution eliminates the interfacial term. One can see that in GaAs/Gao.6Alo.4As 
structures the interfacial correction amounts to less than 10% of the intrasubband Rashba term, and about 20% of 
the intersubband SIA matrix element. 

For a single heterojunction and other structures where the electric field is created by separated charges of impurities 
in depleted doped regions and the 2D electron gas itself, Eq. (fTTjl should be modified by taking into account the electric 
field inhomogeneity. In this case, the electric field £ in the Rashba term should be replaced with its average value 
(integrated with the squared wave function of the corresponding subband). The field inhomogeneity will also contribute 
to the inter-subband term. The correponding corrections can be evaluated using the Poisson equation. For example, 

-t-CXD Z 

the leading correction to the first term in Eq. (fTT|) is: A^^^^ = {ay — ax) {^imge/e) J ^'m(z)5'„(z) f 5'o(-z') dz' dz, 

— oo — oo 

where m ^ n, rig is the sheet concentration of electrons, e is the absolute value of the electron charge, e is the dielectric 
constant, and ^'q is the squared ground-state wave function (here we assume that most electrons are in the lowest 
size-quantization subband; otherwise summation over all occupied subbands should be performed). 
BIA (Dresselhaus) terms in bulk zinc-blende semiconductors have the form Hbia = h**""^ (k) • S, with 

= ah^ (m, ~' K {kl - kj) , (14) 

where irie is the effective mass of the electron. Eg is the band gap, and k^, ky, kz are components of the wave vector 
along the cubic axes [100], [010], and [001] respectively. The y- and z-components of h^""= are obtained from Eq. (fT4|) 
by permutation of indices. 

The expression for the BIA matrix elements between 2D subbands depends on the QW orientation with respect to 
crystal axes. For the most common case of a [100] quantum well 

j^BiA _ ^BiA e _ e ^ , j 



-t-oo 



where 



myj/2mEg 

X J ^^{z)-^^niz)dz{kl-kl)ez, (15) 



+ 00 



A comparison of Eqs. pT|) and (|15p shows that BIA and SIA terms in quantum wells are related to the corresponding 
bulk Hamiltonians in different ways. Intra-subband BIA terms are determined by a unique constant a, the same 
as in bulk BIA terms. This is confirmed experimentally: the values of a measured from spin relaxation in bulk 
GaAs,^^ and from spin-flip Raman scattering in GaAs/AlGaAs quantum wells,'22. are indeed very close (0.07 and 0.065, 
respectively). On the contrary, Rashba terms depend on a specific combination of constants (ay — ax) describing 
different contributions from gradients of electrostatic and crystal potentials, while, for instance, only ay contributes 
to spin-dependent scattering from charged impurities in the bulk. These constants cannot be determined separately, 
only by measuring spin-orbit splitting in size-quantization subbands. These measurements do not, therefore, give 
complete information about the spin-orbit Hamiltonian of the conduction band. 
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Experimental determination of inter-subband spin-flip matrix elements, which would be very helpful in this view, is 
hindered by the fact that they are typically much smaller than the energy of size quantization, and all the observable 
effects of these matrix elements are consequently suppressed. We should like to note that this difficulty can, in 
principle, be overcome by tuning energy levels of quasi-2D electrons with a strong magnetic field. When a magnetic 
field is applied normal to the QW plane, the in-plane motion of electrons is also quantized, and the continuum energy 
spectrum of each 2D subband is transformed into discrete Landau levels (LLs). Landau levels belonging to different 
subbands n and m can be tuned in resonance by choosing the magnetic field so that |i?„ — Em\ = IfkOc, where I is an 
integer, and lOc is the cyclotron frequency. Under these conditions, the energy spectrum and the structure of wave 
functions of the four-level subsystem formed by Zeeman-split spin components of the two LLs from different subbands 
are expected to be strongly affected by the spin-orbit interaction. 

Using the Landau gauge {A^ = —Hy; Ay — ~ 0) and Eqs. (O, (O, (fTTjl . (fT5|) and (fTB|) . we obtain the following 
expressions for matrix elements between electron states differing by subband number, LL number, and spin projection 
on the normal to the structure plane: 

1 



{7J,n,k,,-l/2\Hso\^^,n-l,k,,+^/2)=iAf//-^V2n, (17) 



{rj,n- l,k,,+l/2\Hso\f^,n,k,, ^1/2) ^~iA^'/^V2^, (18) 
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(7?,n,fc,,+l/2|i?so|/x,n-l,/c„-l/2) =A^^^^ — V2n, (19) 

{7^,n~l,h,,~l/2\Hso\f^,n,h,,+l/2) =A^^^^V2^, (20) 

where indices rj and ^ enumerate subbands, n is the LL number, and Ih is the magnetic length. Remarkably, states 
from different subbands and LLs are coupled by either BIA or SIA terms, but not by the two types of spin-orbit 
terms at the same time: Eqs. ([T7)) and (|18p couple states {n — 1, -l-l/2| and (n, — 1/2|, while Eqs. ^9)1 and pp)) couple 
states {n — 1,— 1/2| and (n, -|-1/2|. Therefore, by tuning a specific pair of levels into resonance one can prepare an 
effective two- level system whose spectrum is determined by either the BIA or SIA term. As Eqs. (fT7|) - ([20|l are based 
on general symmetry properties of Rashba and Dresselhaus terms rather than on any specific method of derivation of 
corresponding constants, they can be used for experimental determination of spin-orbit parameters. 

The exact form of this spectrum and its experimental manifestations depend on the electron concentration in the 
quasi-2D system. If the concentration is low and many-body effects are absent, the inter-subband spin-orbit coupling 
results in anticrossing of the two levels. The resulting energy gaps, equal to corresponding inter-subband matrix 
elements (Eqs. ([TTt - ipU)) ') can be measured by high-resolution spectroscopic techniques (optical or spin/cyclotron 
resonance spectroscopy), as has been proposed for one-subband systems in tilted magnetic fields.^^ The admixture 
of wave functions with opposite spin near the anticrossing can, in principle, be detected using optical polarization 
spectroscopy. If all the electron states below the two selected levels arc filled, many-body effects have been shown to 
strongly affect the spectrum of electrons. In particular, they result in the opening of gaps in the energy spectrum 
even in the absence of spin-orbit interaction.— It has been proposed to employ non-equilibrium phonons to probe 
this system of strongly coupled electrons.^"* In this type of experiment, the spin-orbit inter-subband coupling would 
manifest itself by allowing phonon-assisted transitions between levels with different projections of the electron spin, 
which are otherwise forbidden i^i^^. 

In conclusion, we have analyzed the transition from bulk to two-dimensional behavior of the spin-orbit interaction 
in multi-subband quasi-two-dimensional heterostructures. Intra- and inter-subband matrix elements of the spin-orbit 
Hamiltonian are derived, and the role of interfacial effects is estimated. The two types of spin-orbit terms (Rashba and 
Dresselhaus) can be distinguished by measuring the inter-subband coupling of Landau levels in quantizing magnetic 
fields. 

We are grateful to the Editors for inviting us to submit a paper to the special issue devoted to the memory of V.I. 
Perel'. His work and unique personality strongly influenced our research and lit up our lives. 
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